Shear viscous effects on the primordial power spectrum from warm inflation 
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We compute the primordial curvature spectrum generated during warm inflation, including shear 
viscous effects. The primordial spectrum is dominated by the thermal fluctuations of the radiation 
bath, sourced by the dissipative term of the inflaton field. The dissipative coefficient T, computed 
from first principles in the close-to-equilibrium approximation, depends in general on the temper- 
ature T, and this dependence renders the system of the linear fluctuations coupled. Whenever 
the dissipative coefficient is larger than the Hubble expansion rate H, there is a growing mode in 
, the fluctuations before horizon crossing. However, dissipation intrinsically means departures from 

^ ' equilibrium, and therefore the presence of a shear viscous pressure in the radiation fluid. This in 

turn acts as an extra friction term for the radiation fluctuations that tends to damp the growth of 
, the perturbations. Independently of the T functional dependence of the dissipation and the shear 

viscosity, we find that when the shear viscous coefficient C,s is larger than 3pr/-ff at horizon crossing, 
pr being the radiation energy density, the shear damping effect wins and there is no growing mode 
in the spectrum. 
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INTRODUCTION 



Cosmological observations, and in particular Cosmic Microwave Background (CMB) measurements [H-Qj S'l'e 
consistent with a nearly gaussian and practically scale invariant spectrum of primordial perturbations, as pre- 
dicted by the inflationary models [1, Q . An early period of inflation also accounts for the inferred flatness of the 
Universe, and provides a solution to the horizon problem. This makes inflation a robust candidate to account 
I for the early evolution of our universe. Inflationary predictions are characterized by the spectral index of the 
. primordial spectrum, its tensor contribution, and the level of non-gaussianity. Present CMB data however sets 
' at most an upper limit on the level of the tensor contribution and non-gaussianity [l], and it is not yet able to 
\^ , discriminate among the different implementations and models of inflation. This situation is expected to change 
' with the next generation of CMB experiments, like ESA's Planck surveyor satellite [S], which will further improve 
our knowledge of the cosmological parameters. 

Inflation in brief is no more than an early period of accelerated expansion. In the standard picture of inflation, 
denoted as cold inflation, the universe rapidly supercools, and inflation should be followed by a reheating period, 
during which the inflationary vacuum energy is converted into radiation. For reheating to take place, the inflaton 
^ field has to couple to other degree of freedom, such that it can decay into light, relativistic degrees of freedom 
. that thermalize and provide the radiation bath Q. During cold inflation one assumes that those couplings 
play no role during the accelerated expansion. The alternative, called warm inflation [sl, Q (for earlier related 
work see [l^l), assumes on the contrary that those coupling can lead to non- negligible dissipative effects, and 
radiation production can occur concurrently with the inflationary expansion. Both background evolution and 
inflaton fluctuations are modified with the inclusion of an extra friction term T0 accounting for the transfer 
of energy between the inflation and the radiation. The dynamics of the fluctuations are now governed by a 
Langevin equation including a noise force term from the influence of the radiation fluctuations into the inflaton 
field [1, [TT| - [l3 | . Thermal fluctuations in the radiation are transfered to the inflaton and become the main source 
of primordial fluctuations ^, JJ., i5,] . 
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The dissipative coefficient can be computed from first principles i n q uantum field theory, within an adiabatic 
approximation. The two-stage interaction configuration proposed in [ISj has been shown to lead to a large enough 
dissipative coefficient while keeping the corrections to the infiationary potential under control and allowing a 
period of slow- roll inflation |16l4l8| . For the two-stage mechanism, the inflaton field couples to a heavy catalyst 
field, and the latter in turn couples to light degrees of freedom. During the motion of the background inflaton, it 
excites the catalyst fields which then decay into light fields [l^ . Using this mechanism, the first calculation of the 
dissipative coefficient within the close-to-equilibrium approximation was done in [2^, leading to a temperature 
dependent dissipative coefficient. In the low-temperature regime, when the mass of the heavy catalyst field is 
much larger than T, one has T oc T"^, while in the high-temperature regime, when T is above the heavy catalyst 
field mass, dependingon type of interaction the dissipative coefficient becomes linear with T [2l[ or goes goes 
like the inverse of T [22[ . All high-T models suffer in general from very large thermal corrections which spoil the 
flatness of the potential [1^ [2J|, with only a few exceptions [H, US]- However, viable models of warm inflation 
have been studied in the low-T regime |25l - l27| . 




The temperature dependence of the dissipative coefficient induces the coupling of the field and radiation 
fluctuations as shown in [2^ . Previous studies of the primordial spectrum of perturbations in warm inflation [l^ 
did take into account the influence of the thermal fluctuations on the field through the noise term, but not the 
coupling through the dissipative term itself. In (28( 1 ^^.s shown that for positive power of T in T, and when 
T dominates over the Hubble expansion rate, this coupling induces a growing mode in the fluctuations before 
horizon crossing, enhancing by several order of magnitude the amplitude of the primordial perturbations with 
respect to previous calculations. In the calculations of the primordial spectrum typically the radiation bath is 
treated as a perfect fluid, with an equilibrium pressure Pr — Pr/'^, where pr is the radiation energy density, an 
approximation valid in the close-to-equilibrium regime required for the calculation of the dissipative coefficient to 
hold. However, even in that regime, the radiation bath is expected to depart from an ideal fluid as a consequence 
of the constant production of radiation particles from the background held dissipation. Imperfect fluids have 
dissipative effects that can be parameterized in terms of shear and bulk viscosity coefficients, and a heat flow 
coefficient. Heat flow happens as a consequence of changes of conserved charges other than the stress-energy 
tensor, but we do not consider this possibility here and focus on the effects of the temperature. Bulk viscous 
effects, which can be interpreted as a consequence of the decay of particles within the fluid, have been considered 
for warm inflation in (29j , where they studied either a constant bulk viscous pressure or one proportional to the 
radiation energy density. The bulk pressure appears at both the background and the perturbation level, and being 
a negative pressure, it will favor warm inflation. For the amplitude of the spectrum, for the phenomenological 
model considered in [2^ they found that it could induce a variation in the amplitude of the order of 4%. On the 
other hand, the shear viscosity is related to changes in momentum of the particles of the fluid, and appears only 
at the level of the perturbations. Shear and bulk viscosity coefficients due to light field have been extensively 
compute in the literature [30l[3l|. leading to power- law dependences on T for these coefficients. In addition, the 
bulk viscous coefficient typically ends being much smaller than the shear viscosity. Therefore, we will concentrate 
on the effect of the shear viscosity on the spectrum, and do not consider those of the bulk viscosity. The shear 
viscosity, being related to dissipation, appears in the radiation fluid equation as a friction term which tends to 
damp the growth of the fluctuations [28|. Eventually, the damping effect will dominate over the enhancement 
induced by the dissipative source term. The aim will be therefore to quantify, in a model independent way, this 
effect on the spectrum, and when it will render the system fleld-radiation effectively decoupled. 

The outline of the paper is as follow. In section II we review the basic of warm inflation at the background 
level. In section HI we present the equations for the fluctuations of the coupled fleld-radiation system, when the 
radiation is taken as an imperfect fluid. The numerical solutions for the fluctuations are presented in section 
IV, with and without the shear viscosity. We also study the inflationary model dependence of the results on the 
spectrum by considering two generic model of inflation: a chaotic model with general power p, and a standard 
quadratic hybrid model. In section V we summarize our flndings. 



In any particle physics realization of the inflationary framework, the inflaton is not an isolated part of the 
model but it interacts with other flelds. These interactions may lead to the dissipation of the inflaton energy 
into other degrees of freedom, such that a small percent of the inflaton vacuum energy is transferred into other 
kinds of energy. In the two-stage mechanism for warm inflation, dissipation leads to particle production of light 
degrees of freedom. When those relativistic particles thermalize fast enough, say in less than a Hubble time in 
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an expanding universe, we can model their contribution as that of radiation: 



(2.1) 

where T is the temperature of the thermal bath, and the effective number of light degrees of freedom^ . 
The dissipative term appears as an extra friction term in the evolution equation for the inflaton field 0, 

0+(3i? + T)0 + 1/0 = 0, (2.2) 

T being the dissipative coefficient, H = a/a is the Hubble rate of expansion, and a the scale factor of the 
Friedmann-Robertson- Walker background metric: 

ds'^ = ~df + a{tf5,jdx'dx^ . (2.3) 

Eq. (j2.2p is equivalent to the evolution equation for the inflaton energy density p^: 

p^ + 3H{p^+p^)^^T{p^+p^), (2.4) 

with pressure = 0^/2 — y ((/>), and pc/, + Pcj, — <p^- Energy conservation then demands that the energy lost of 
the inflaton field must be gained by the radiation fluid pr , with the RHS of Eq. (12. 4p acting as the source term: 



Pr + 'SH{pr+Pr) ^T{p^+p^). (2.5) 

In warm inflation, radiation is not redshifted away during inflation, because it is continuously fed by the inflaton 
through the dissipation [9]. Inflation happens when p/j <C p^, but even if small when compared to the inflaton 

energy density it can be larger than the expansion rate with p^J^'^ > H . Assuming thermalization, this translates 

roughly into T > H. Otherwise, when T < H (or similarly when p^^* < H), one just recovers the standard cold 
inflation scenario, where dissipation can be neglected. 

During warm inflation the motion of the inflaton field has to be overdamped in order to have the accelerated 
expansion, but now this can be achieved due to the extra friction term T instead of that of the Hubble rate. And 
once </>, H , and also T, are in this slow-roll regime, the same will happen with the radiation energy density, the 
source term now compensating for the Hubble dilution. In the slow-roll regime, the equations of motion reduce 
to: 

3H{1 + Q)<P~~V^, (2.6) 
4pR c 3Q02 ^ (2.7) 

where we have introduced the dissipative ratio Q = T/{3H). Notice that Q is not necessarily constant. The 
coefficient T will depend on (j) and T, and therefore depending on the model the ratio Q may increase or decrease 
during inflation [25]. The slow-roll conditions are given by [18]: 

S = ^<1. (2.11) 

y<h 



^ If not otherwise specified, we will take = 228.75, the effective no. of degrees of freedom for the Minimal Supersymmetric 
Standard Model, when presenting numerical results. 
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where the slow-roU parameter /3x takes into account the variation of T with respect to (j), and the last condition 
ensures that thermal corrections to the inflation potential are negligible. Similarly, taking also into account the 
dependence on T of T, one has: 



dlnT 



dlnT 



<4, 



(2.12) 



which reflects the fact that radiation has to be produced at a rate larger than the redshift due to the expansion of 
the universe. The slow-roll regime ends when any of the above conditions Eqs. (|2.8I) - (|2.11I) is no longer satisfied, 
such that either the motion is no longer overdamped and slow-roll ends, or the radiation becomes comparable to 
the inflaton energy density. Either way, inflation will end shortly afterwards. 

For warm inflation the first requirement is to have T > H, but the ratio Q can be larger or smaller than unity. 
In the former case we are in the strong dissipative regime (SDR), whereas the latter is called weak dissipative 
regime (WDR). In the weak dissipative regime the extra friction added by T is not enough to substantially modify 
the background inflaton evolution, and it will resemble that of cold inflation; still the thermal fluctuations of the 
radiation energy density will modify the field fiuctuations, and affect the primordial spectrum of perturbations. 

The T and </> dependent dissipative coefficient has been computed in [20[, using the near-equilibrium approx- 
imation proposed in [2^ . The specific field theory models considered for the inflaton interactions leading to 
dissipation all follow from the two-stage mechanism [l^, [l^, [s^l • In this mechanism, the inflaton field (j) is cou- 
pled to heavy catalyst fields x, which decay into light fields at. Consistency of the approximations then demands 
the microphysical dynamics determining T to be faster than that of the macroscopic motion of the background 
inflaton and the expansion: 



(2.13) 



where is the decay width of the heavy flelds. In addition, the condition T ^ H allows to neglect the expansion 
of the universe when computing T. In the low T regime, when the mass of the catalyst field is larger than 
T, one has: 



2^3 ^3 

T{(f>,T) fx — oc — 



(2.14) 



while in the high T regime, where the thermal corrections to the catalyst field mass start to be important, 

T(0,T)cxr. (2.15) 

And in the very high T regime, the dissipative coefficient goes like the inverse of T. These are the cases of study 
we are going to consider in the next section when studying the fiuctuations during warm infiation. In general, 
we will parametrize the dissipative coefficient as: 



T = C, 



(2.16) 



with c — TO = 1. We will work with c =3,1,-1, and also c = 0, the case of no T dependence for the dissipative 
coefficient. 



III. FLUCTUATIONS AT LINEAR ORDER: PRIMORDIAL SPECTRUM 



During warm inflation we have a multicomponent fluid, a mixture of a scalar inflaton field $ interacting 
with the radiation fluid. Both components exchange energy and momentum through the dissipative term T. 
Dissipative effects also imply small departures from equilibrium, and that the radiation fluid will not behave 
exactly like a perfect fluid during inflation. In relativistic theory, these effect can be parametrized in terms of a 
shear viscous tensor TTgh, an energy flux vector qa and a bulk viscous pressure tt^, in the stress-energy tensor for 
the radiation fluid (ssl , 13^ , 

T^l^ = (Pr +Pr+ ^bW:K'''> + iPr + ^b)9ab + Q^I^ u['^ + q^^^^^^ + TT.fe , (3.1) 
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_ . . _ . . (r) 

where pr is the energy density, Pr the adiabatic pressure, Ua the four velocity of the radiation fluid, gab the 
four-dimensional metric, and u''a\°''' = = gabT^""^ i u^a^'q"' = 0. There would be heat flow for example in the 
presence of conserved charges in the system other than the stress-energy tensor, but we do not consider such 
possibility in this work, and then set = 0. The shear viscous tensor vanishes in an homogeneous and isotropic 
background geometry, but at linear order it is given by [33 |: 

T^ab — -"^CsC^ab , (3-2) 

where is the shear viscosity coefficient and Gab the shear of the radiation fluid: 

Oab = V(QUfc) + Ui^aU'VcUb) - ^^VUc , (3.3) 

Va being the covariant derivative of the metric gab- The bulk viscous pressure can be seen as a non-adiabatic 
pressure contribution, already present at the background level. Nevertheless, the contribution from the light 
fields is expected to be small with respect to Pr- Thus, we will also set tt;, = and focus on studying the 
consequences of the shear viscosity on the primordial perturbations during warm inflation. 

In order to study the system of perturbations at linear order, field, radiation energy density and radiation 
pressure are expanded around their background values in a Friedman-Robertson- Walker metric: 

<^{x,t) = (t,{t)^5(t){t,x), (3.4) 

Pr{x,t) ^ Pr{t) + 5pr{t,x) , (3.5) 
Prix,t) = Pr{t) + Spr(t,x) , (3.6) 

and similarly for the dissipative coefficient: T{x,t) — T(t) + 6T{t,x). The perturbed FRW metric, including 
only scalar perturbations, is given by^: 

ds'^ = -(1 + 2a)dt^ - 2ad,fidx'dt + a^[5^j{l + 2ip) + 2didj-/]dx'dx' , (3.7) 

where tp is the intrinsic curvature of a constant time hypersurface. For later use, we introduce the combinations: 

X = a{l3 + aj), (3.8) 

K = 3{Ha^^) + dkd''x, (3.9) 

where x is the shear and —k the perturbed expansion scalar of the comoving frame. 
The evolution equations follow from the conservation of the energy-momentum tensor: 

V'Ti-^^Qi"), E^?i"'-0, (3.10) 



where Qb is the four-vector source term accounting for the exchange of energy and momentum 

- ( 

is now the four-velocity of inflaton fluid: 



(0) ^ gW ^ Y^.;Vaa>V6$ , (3.11) 



(3.12) 



and then: 

QI'^ = T(p^+p^)i/2v,$. (3.13) 

The projection of the four- vector source term along the direction of the fluid gives the energy density source 
term, 

Q^''^ = -ulQi^^ , (3.14) 



^ Latin indexes i, j, k, . . . are used for the spatial components, and either a, b, c, . . . or Greek letters for space-time indexes. 
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which at linear order is given by: 

= + <5Q, , (3.15) 

Qr = T4>^ , (3.16) 

6Qr ^ 6T<j}^ + 2T^6<})~2aT^^ . (3.17) 

The momentum source term Ja is orthogonal to the fluid velocity: 

Q(f = Q^^^ui"^) + , = , (3.18) 

and vanishes in the FRW geometry; at linear order it reads: 

jf ^ = d,3r , (3.19) 
Jr = -T^5(j). (3.20) 

To complete the specification of the source, we need ST, which for a general temperature T and field (p dependent 
dissipative coefficient, T = Cci,T'^/(j>™, with c — m = 1, is given by: 

^T = T(cf -^1). (3.21) 

Although dissipation implies departures from thermal equilibrium in the radiation fluid, the system has to be 
close-to-equilibrium for the calculation of the dissipative coefficient to hold, therefore pr ~ Pr/3, Pr (x T"* and 

4-^. (3.22) 

T Pr 

Finally, the evolution equations for the radiation fluctuations are obtained expanding at linear order Eq. p.lOp 
[ssl - lssj . Working in momentum space, defining the Fourier transform with respect to the comoving coordinates, 
the equation of motion for the fluctuations with comoving wavenumber k are given by^: 

6pr + 3H{dpr + Spr) — -3{pr + Pr)ip -i j + {pr + Pr)x] + 5Qr + QrOi , (3.23) 



+ iH^r = -{pr + Pr)a - Spr + — CT^ + Jr , (3.24) 



2fc2 
3a 

where a "dot" denotes the derivative with respect to the metric time "t" , is the radiation momentum pcrtur- 



bation, T*^ —dj'i'r/a, and cr^ the shear viscous pressure at linear order: 



ar :^ -2Cs + x) ■ (3.25) 

\Pr+Pr J 

On the other hand, field fluctuations are described by a stochastic system whose evolution is determined by a 
Langevin equation [2^, : 

S4> + (SH + T)4 + (^^ + 1^00) H = [2(T + H)Tf^ a'^/^^k - ST^ 

+ (pin + a) + (20 + iH^)a - T{5p - ap) , (3.26) 
where the stochastic source can be approximated by a localized gaussian distribution with correlation function: 

{^{t,x)^(t',x'))^5{t-t')5'^^\x-x'). (3.27) 



^ For simplicity, we keep the same notation for the fluctuations (5/{x, t) and their Fourier transform 5/(k, t). 
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So far, the equations for the perturbations at linear order are written in a "gauge ready" form, without 
specifying any particular gauge, but the equations can also be written in terms of gauge invariant (GI) variables. 
For any scalar quantity /, at linear order one can define a gauge invariant perturbation (ssl |36|: 



Sf""' = 6f-I^^, (3.28) 



while similarly the gauge invariant momentum perturbation reads: 



H 

and for the metric perturbations one has the gauge invariant combinations: 



M,Gi = ^_P±P^^ (3.29) 



(3.30) 



-Hi 

<i> = ip-Hx. (3.31) 



The evolution equations then read: 

a 



T0yi + M + 2(0 + 3i/0)yi- ^</>— , (3.32) 

H 



h-2 

<5/j^^ + 3i7(l + wr)5p'^^ = + <50^^ + PrA , (3.33) 

/ i.2 _ \ oiL2 _ ^ 

*f + 3i/ 1 + -^Qs = -WrSp'^' - TcbScj)''' - {pr+Pr)A - —{p^+p,)^- , (3.34) 

where in Eq. p.34p we have defined: 

y Pr+Pr 

Finally, from the Einstein equations at linear order, the gauge invariant metric perturbations are given by [35l.[36|: 

A = -^:R, (3.36) 

fc2 3Sp^^ , , 

$ = 3^1 + --^, (3.37) 



a2i?2 2 pT 

where !K is the total comoving curvature perturbation, 

= LP "^T 

PT +PT 

H 



. (3.38) 

PT +PT 

For a multicomponent fluid, the total momentum perturbation is given by the sum of the individual components, 
and IK can be written as the weighted sum of the individual contributions: 

^= y. 7^3?„, (3.39) 
3ia = -#-*^^ (3.40) 



where we have defined = Pa + Pa- In particular, for a scalar field = —(f)5(fP^^ — 0^, and 

H 



= ^6<j)G' , (3.41) 
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for which the power spectrum would be the form for single field cold inflation. For warm inflation, we shall use 
the total comoving curvature perturbation to evaluate the primordial spectrum, 

P^{k)^^{\-Rk?), (3.42) 

where "(•••)" means average over different realizations of the noise term in Eq. p.26p . In all numerical results 
shown in this work we have performed averages over 1000 runs. This was found to be more than enough to get 
convergent numerical results." . 

The largest observable scale in the CMB corresponds to a comoving scale k crossing the horizon e-folds 
before the end of inflation, denoted by A: = a^H^. The value of Ng can vary in general from 40 to 70 depending on 
the inflationary model and details on the subsequent reheating process Although we shall consider different 
inflationary potentials, we will not consider the details of reheating period, and simply fix the horizon crossing at 
Ng ~ 60. As we will see, soon after horizon crossing, the amplitude of the individual curvature perturbations 
and 5ir freezes out, and so does that of the total curvature. This allows to compute the primordial spectrum by 
evaluating Eq. (j3.42p at horizon crossing, mainly when getting analytical approximations. Nevertheless, when 
showing numerical results we shall integrate Eqs. p.32p - p.34p and evaluate the amplitude of the spectrum say 
10 e-folds after horizon crossing. 



IV. EQUATIONS IN THE ZERO-ORDER SLOW-ROLL APPROXIMATION AND BEYOND 



In order to gain some insight on the evolution of the perturbations, we follow 28] and first study the equations 
for the perturbations at zero-order in the slow-roll parameters. That is, expanding the background variables 
around their slow-roll values and neglecting all terms in the equations proportional to slow-roll parameters. This 
eliminates in the equations for the fluctuations the dependence on the inflationary potential, and the details 
on the evolution of the dissipative coefficient. For example, the metric perturbation A given in Eq. (I3.36P is 
proportional to the slow-roll coefficient e, Eq. (|2.8I) . e = —H/H^, and can be neglected; and the same for the 
terms proportional to H(f), i.e., those proportional to <&, and (j)/H(j). Defining the dimcnsionless variables: 



Vk 



Wk 



Uk 



k^/^d<j>' 



GI 



[2{T + H)T]^''^ ' 



[2(T-i-i/)r]^/^0 ' 

and using the slow- roll background equations p.6p . p.7p . we have the system: 



yk + iH{l + Q)yk + H^ 



3r?(l + Q) - 3mQ 



Vk 



3/2 



^fc - iQcH Wk , 



Wk + H{A - c)wk = ^^^"fe + ^yfc , 



Vk 



ilk + SH{1 + z^Cs)uk = ~3QH (- 
where z = k/(aH). Combining Eqs. (|4.5p and (|4.6p into a second order differential equation, we have: 



' Wk 

3 



Wk + H{9 - c + 3z^Cs)wk + 



20-5c + 6Qc+ — + 32^(4 - c)Cs 
o 



Wfc = 2 - 
a 



3/2 



Hi4.-6Q + 6z^Cs)yk-H^ 



3z^ + 6T]il + Q) - 6Qm-^ 



Vk ■ 



(4.1) 
(4.2) 
(4.3) 



(4.4) 

(4.5) 
(4.6) 



(4.7) 
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(aH/k) 

FIG. 1: Evolution of the power spectrum of the field {yi), radiation energy density (wl), and radiation momentum 
(^fe)/Q'^i for Q ~ 100, and wavenumber k — 10*ff. The vertical thin dotted line sets the value of the freeze out scale kp 
in warm inflation. The results are shown for different power dependence on T of the dissipative coefficient: c = 3 (solid 
lines), c = I (dashed lines), c — ~1 (dash-dotted lines), and c — (dotted lines). 




(aH/k) 

FIG. 2: Evolution of the total curvature perturbation spectrum P^^^ (black lines), the radiation P^^^ (red lines), and the 
field P^^^ (green lines) curvature perturbation spectrum. The results are shown for different power dependence on T of 
the dissipative coefficient: c — 3 (solid lines), c = I (dashed lines), c = — f (dash-dotted lines), and c = (dotted lines). 

Before including shear effects, we study the evolution of the perturbations setting — 0. The evolution of 
the power spectrum of the field (j/^), the radiation energy density (w^), and the radiation momentum (m^), are 
shown in Fig. [T]as a function of z^^ — aH/k. We have taken Q = 100, and started the integration at Zi = 10^. 
As mentioned above, quantities as (y^) denotes the average over 1000 realizations of the gaussian noise term, 
and by y^, w^, u\ we mean the modulus of the complex variable. We have set the initial conditions for field 
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fluctuations in the vacuum, while Wk and u^, initially vanish, for simplicity. Starting the evolution early enough 
before horizon crossing, the system is always controlled by the stochastic source term, and the dependence on 
the initial conditions is quickly erased. 

We have considered different powers of T for the dissipative coefhcient T, c = 3, 1,-1, and included the case 
of a constant or field dependent T { c — 0, dotted lines) for comparison. The radiation fluctuation Wk acts 
as a source term for the field, but at early times ^ 1, for subhorizon perturbations, the field evolution is 
dominated by the stochastic source term, and both radiation and field fluctuations evolve like in the case c = 0. 
In the latter case, the freeze out of the perturbations takes place before horizon crossing, due to the extra friction 
term in Eq. (|4.4I) . at around kp/{aH) ~ 3y^Q/2 (vertical thin dotted line) [12, 15, 28], and soon after field and 
radiation spectrum level off. The field spectrum for a T independent T can be computed analytically and is 
given by [l2| : 



where the subindex "0" denotes the value for c = 0. On the other hand, when c > 0, field and radiation 
fluctuations get effectively coupled before freeze out at around « 18Qc, and both start growing at similar 
rates. Numerically, we get: 



The field spectrum when the fluctuations are still subhorizon and z > z^, which is independent of the radiation 
fluctuation, can be found in .28.] : 



and therefore, at horizon crossing: 



ivl)^.^-, (4.10) 

Zc 



(y^).oczf-iocQ(^-i)/^ (4.11) 



When c < 0, the effect is the opposite, and effectively the freeze out is delayed until practically horizon crossing, 
which makes the amplitude of the field spectrum smaller than in the c = case. 

Fig. [2] shows the evolution of the comoving curvature power spectrum, Pj^ , given by the sum of the radiation 
and the field contributions as in Eq. p.39p . for different values of c. Also shown are the power spectra of the 
radiation, , and the field Py/^ ■ After horizon crossing they all converge to the same amplitude. During 
slow-roll one has that hr = Pr + Pr — Qh^, and from Eq. (|4.6p when z <C 1 the radiation momentum becomes 
proportional to the field fiuctuation; 



^r^Q^S(j), (4.12) 



and therefore: 



Pjl,. = YP^r- ■ (4.13) 

Owing to the fact that <^ hr ~ /it, the main contribution to the primordial spectrum in Eq. p. 421) comes 
from the radiation, before and after horizon crossing. The primordial spectrum is always dominated by the 
thermal ffuctuations. But after horizon crossing one simply has: 

^ Pa?,, ^ Pk^ . (4.14) 

Therefore, the amplitude of the primordial spectrum can be written as usual in terms of that of the inflaton field: 

fHV (H + T)T , 

^2 (yl)*^ (4-15) 

evaluated at horizon crossing. 
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FIG. 3: The spectrum of the field at horizon crossing as a function of the dissipative parameter Q, at zero order in 

the slow-roll parameters, for different values of c, and no shear = 0. Solid lines are obtained integrating Eqs. (|4.4|l - (|4.6|l . 
while dashed lines were obtained with the approximation used in ^2^. 



c 


a 




A, 


Be 


3 


7.5 


7.0 


1.9 X 10"** 


3.4 X 10"'' 


1 


2.5 


2.0 


2.8 X 10"^ 


6.8 X 10"^ 


-1 


0.2 


1.4 


0.78 


0.088 



TABLE I: Coefiicients for the numerical fit of the spectrum, Eq. (|4.16p and Eq. (|4.17p . 



In Fig. [3] we have compared the power spectrum of the field, {yf} at horizon crossing, as a function of the 
dissipative parameter Q for different values of c. The equation for the fluctuations has been integrated keeping the 
background values constant (solid lines), Eqs. ()4.4|) - (|4.6p . The larger the power c > 0, the larger the enhancement 
with Q, as the fluctuations get coupled earlier. For c < 0, as mentioned before, we have the opposite effect, and 
the spectrum diminished with respect to the case c = 0. For positive c, the curves can be fitted by a function: 

{yl).\c>o ^ ivlMAcQ" + B,QP) . (4.16) 
but when c = — 1, we have found that the curve can be best fitted by: 

The coefficients are given in Table H For c = 3, 1, the approximation works well for Q > 50, while for c = — 1 it 
is valid for any Q. 

We have also included in Fig. |3]the spectrum of the field obtained with the approximation used in [28l[ (d ashed 
Hnes) for comparison. We have confirmed the main results about the power spectrum obtained in |28 |: i.e., 
that for c > 0, the amplitude of the primordial spectrum in warm inflation is enhanced by a factor ~ 
However, while in [28| they found a = 3c, we have a smaller power a — 5c/2, which for c ~ 3 can mean a 
difference of two or 3 orders of magnitude in the amplitude of the spectrum for Q ~ 50 — 100. The differences 
can be traced back to how the radiation source term is treated. In ; 2^ . it was approximated by: 

SQr^Qr^C^cQr^, (4.18) 

while we have kept the dependence on the field: 

<5Q.-g.fc^ + 2^) . (4.19) 
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FIG. 4: The spectrum of the field {yf.) at horizon crossing as a function of the dissipative parameter Q,, for different 
inflationary models, c = 3 and no shear C,a ~ 0. The solid line is the result at zero order in slow-roll; dashed lines are for 
a chaotic model with p =6,4,2, from top to bottom; the dash-dotted line is a quadratic hybrid model. 



With no shear viscosity included, the second order differential equation Eq. ()4.7p reads: 



H^GQcwk + H{A - 6Q)yk - 



3/2 



22^ + 677(1 



Q) - 6Q™-^ 



(4.20) 



where we have written the equation such that on the RHS we have the terms induced by the field dependence 
in Eq. (|4.19p . Thus, setting the RHS to zero one recovers the equation derived with the source term as given in 
Eq. (|4.18|) . While the term proportional to the field perturbation i/k acts as a source term on the radiation that 
tends to enhance the fiuctuation, the extra terms GH^QcWk and iJ(4 — 6Q)yk have the opposite effect, i.e., that 
of damping the growth. Although these contributions are not enough to avoid the growth of the fluctuations, 
they have a sizable effect on their power-law behavior, mainly when c = 3, as seen in Fig. [3] When c = — 1 
the radiation fluctuations do not grow before horizon crossing, so that the effect of the fleld dependent terms in 
Eq. (|4.19p is negligible. 

Neglecting the evolution of the background variables and working at zero order in the slow roll parameters, is 
a good approximation when Q is large enough, and the background parameters hardly vary during the last 60 
e-folds of inflation. Indeed the spectrum depends mainly on the value of the parameters in a smaller interval, 
5-6 e-folds around horizon crossing, where one may expect the approximation of keeping them constant to be 
a fairly good one. Still, this is a model dependent question. This can be seen in Fig. 21 where we show the 
field spectrum for some generic inflationary models, and c ~ 3, Q = 0. The value of (?/^.)* has been obtained 
integrating Eqs. p.26p - p.34p . evaluating the amplitude of the comoving curvature spectrum at Ne = 20 efolds 
after horizon crossing, and using Eq. (|4.15|) . We have considered the inflationary models: 



V 



V = Vo 



p \mp 



2 \mp 



(chaotic) . 
(hybrid) . 



(4.21) 
(4.22) 



For the chaotic model we have run different powers p = 6, 4, 2, and set Vq 



10- 



while for the hybrid 



Vb = 10 ^nip and 77^ — 3. For each model, the initial value of the inflation field is chosen such that we can have 
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c 


As 


Bs 


A, 


3 


6.35 


3.4 


1.36 


1 


1.9 


1.2 


1.33 


-1 


-0.95 


-0.7 


-0.66 



TABLE II: Coefficients for the numerical fit of the spectrum of the field with the shear viscosity, Eq. (|4.23p . 



A'e — 64, and from the background slow-roll equations one derives the initial values of the field derivative, p,- and 
Q. We have chosen k = lOOi/i, Hi being the initial value of the Hubble parameter. Therefore, horizon crossing 
k = at,H^ takes places at around 60 e- folds before the end of infiation. 

In all the examples considered above the dissipative coefficient increases during infiation 25] , and the larger the 
power in the potential, the slower the evolution of the background values. In the plot, the solid line is the result 
obtained with constant background variables. For a quartic chaotic model or larger power, the approximation at 
zero order in the slow-roll works fine, while for a quadratic chaotic it tends to overestimate the spectrum by at 
least an order of magnitude for Q > 50. For the hybrid model, the model dependence shows up when Q < 100. 

Shear effects will further damp the growth of the fiuctuations. In Eq. (|4.7I) the shear acts as an extra friction 
when the fluctuations are still subhorizon, suppressing the amplitude of the radiation fluctuation before the 
radiation-field system becomes effectively coupled. Whenever the shear is large enough, this suppression indeed 
can prevent altogether the growth of the field perturbations, as the amplitude of the radiation fluctuation is not 
enough to affect that of the field before horizon crossing. This will happen when Cs = C,sH/ {2>pr) ^ 1 at horizon 
crossing. During warm infiation, we have the catalyst field coupled to the inflaton field, and to the light degrees of 
freedom giving rise to the thermal bath. The shear viscosity for light fields (light with respect to the temperature 
T of the thermal bath) typically behaves as oc [Hillol, although, depending on the pattern of interactions, 
other powers could be possible and cannot be excluded. Nevertheless, the damping is fully controlled by the 
value of the dimensionless parameter Qs at horizon crossing, and is independent of the functional form of the 
shear with the temperature, as can be seen in Fig. [S] We have integrated the full EOM without approximations 
Eqs. p.26p - p.34[) . for a quartic chaotic model, Eq. (|4.2ip with p = 4, and set the initial background values such 
that Q* — 40. We show in the plot the dependence of the field spectrum with C, at horizon crossing, for different 
values of c, and two examples of the shear viscosity: one proportional to (solid lines), and another linear in 
T (dashed lines), however the curves lay on top of each other. We have checked that this is independent of the 
infiationary model considered. We have normalized the field spectrum with the value obtained when c — Q. As 
the shear increases, it does damp the radiation fiuctuation enough for the field spectrum to approach the c = 
value. Asymptotically, when ^ 1, for a linear dissipative coefficient with T one practically recovers the c — Q 
case, for a cubic one the field spectrum is ~ 2(y^)o/5, while the inverse T case is shghtly above ~ 5{yl)o/2. In 
all cases, the field spectrum is well fitted by a function: 

logio TTT -As-Bs[l + tanh (logio C + A.)] , (4.23) 

which interpolates between the result with no shear ~ 10"*= and the c = case, modulo a normalization constant 
~ 10'^=~^^». As an example, the coefficients Ag, Bg and for the potential shown in Fig. [S]with = 40 are 
given in Table [Til 

Finally, combining Eq. (|4.23p with Eq. (j4.16p . the field spectrum when c > reads: 

{yl),^{yl)oFQ[Qf^^^'\ (4.24) 



where: 

Fq[Q] ^ {A,Q" + B,Q^) , (4.25) 

- ^ [1 - tanh (logio Cs + A,)] . (4.26) 

Therefore, in the strong dissipative regime when Q > 1, the primordial spectrum of the curvature perturbation 
can be written as that obtained for a T independent dissipative coefficient, times an enhancement function -Fq[(3] 
depending on the dissipative ratio Q, but controlled by a function depending on the shear i^s [Cs] ■ 



^ " ' / 4^2 



3/2 



(1 + Q) 

1 + 30 



^) xFQ[g]^=[f=l, (4.27) 
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FIG. 5: Field spectrum normalized by the value with c = 0, as a function of the shear parameter ^s, for different values 
of c, and Q = 40. From top to bottom, c = 3, 1, —1. For each curve, we have also considered two different T dependence 
on the shear viscosity, as indicated in the legend, but both gives the same field spectrum. 

and whenever > 1 one recovers the amplitude of the primordial spectrum obtained when c = 0. The latter is of 
a magnitude compatible with the observational value, for model parameters values common in inflationary model 
building (25l. [26|. For example, without the enhancement, a quartic warm chaotic model gives rise to the right 
level of perturbations with a coupling constant A ~ 10~^^ — 10"^'*. From the point of view of model building, 
the enhancement produced by the backreaction of the radiation fluctuations onto the fields, if not avoided by 
shear effects, can be compensated by lowering the height of the potential, i.e., by lowering couplings and masses. 
However, it will also have an impact on the prediction for the spectral index: for models with Q increasing at the 
time of horizon crossing it will render the spectrum too blue-tilted, and the other way round, for Q decreasing 
the spectrum may become too red-tilted. Shear viscosity damps the growth of the fluctuations, and therefore 
will also affect the spectral index. 

Even though the results we have obtained are fairly model independent and shown not to depend on the specific 
dependence on the temperature of shear viscosity, it is useful to estimate the magnitude of the ratio C.sH/{3pr) 
for typical warm inflation models. In kinetic theory, the shear viscosity for relativistic fluids can be expressed 
parametrically as proportional to the mean free path of quasiparticles in the fluid. Considering as an example a 
radiation fluid made of relativistic scalar particles a, with mass rria/T <C 1 and self-interaction potential Ao-cr'^/4!, 
the mean free path of quasiparticles is determined by the inverse of the thermal width, which is O(A^), and the 
computed value for the shear viscosity is [13] (s — 3 x lO^T'^/A^. The condition (sH/{3pr) > 1 can then be 
expressed, for example, as a condition on the magnitude of the radiation bath self-interaction, Ao- ^ 55y/H/T. 
Since warm inflation requires T > H and in general we typically work with values T ^ H, we find that weakly 
interacting radiation fluids can easily have shear viscosities of sufficient magnitude to counterbalance and suppress 
completely the growth of fluctuations caused by the coupling of the inflaton's fluctuations with those from the 
radiation. 

In supersymmetric warm inflation models however the radiation bath self coupling is the same as the coupling 
to the catalyst fleld, which enters in the calculation of the dissipative coefficient. In general, large multiplicities 
of the fields (catalyst x and radiation a) are required in order to have enough dissipation. Considering a model 
with [M^, [Mo- copies of the complex fields, with common coupling h among the x's and cr's, the self-interaction 
potential is given by /i^ J2i In the low-T regime, the dissipative and shear coefficients can be written as 
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FIG. 6: Field spectrum as a function of Q*, for a quartic chaotic model, with a cubic dissipative coefficient (Eq. I4.28P 
and (s given by Eq. (|4.29|) . We have taken ^f^ = 10^ and h = \^ (solid line), h = 1.67 (dashed hne) and h = 1.12 
(dot-dashed line). We include for comparison the result with c = 0. 



T~0.1/i%:K^ J, (4.28) 
C^l27N^{l + 0.3h)^, (4.29) 

where we have included the next-to- leading order correction in the shear viscosity ^4^. The condition Cs > 1 
then reads: 

1^ J-f J-f 

/i'<128— -^ 69-, (4.30) 

g* T T 

where g^. ~ ISNct/S. Although in a bath of weakly interacting radiation particles, ft. ^ 1, the shear viscosity 
would be rather large, it may not give rise to enough dissipation to sustain a period of warm inflation, unless 
the weakness of the coupling is compensated by having a large multiplicity for the fields. This therefore becomes 
a model dependent question depending on the parameters Jvf^, J^a- and the coupling h. In Fig. [5] we show an 
example for a quartic chaotic model, where we compare the value of the field spectrum (y^)* as a function of 
for different values of the Yukawa coupling h. We have kept the value 3Nfo- = 10'^ fixed, and vary the value of Ji^ 
to get different values of Q,. From top to bottom the value of h decreases, starting with the largest possible one 
h — (solid line) for which there is a negligible shear effect. In this case, a value Q* ~ 10 requires INf^ = 35, 
while for Q* — 100 we need Ji^ ~ 110. As the value of h decreases, the multiplicity of the field to get the same 
value of increases by a factor (47r//i^)^. For h = 1.12 and — 100, we would need 3Nf^ ~ 10^. As the value 
of Q* increases, so it does the ratio T/H, and therefore the parameter (s decreases along the curves. Although 
in this example viscous effects are not enough to completely avoid the growth of the perturbations, they bring it 
down to {yD* oc Q^-^ for h — 1,12 and {yD* oc Q^-^ for h — 1.67, instead of {yD* oc . 



V. CONCLUSIONS 



Density perturbations in warm inflation are seeded by thermal fluctuations of the inflaton. In warm inflation 
the inflaton decays into radiation through a dissipation term in the inflaton's equation of motion and that 
originates from the microphysical interactions of the inflaton fleld with other degrees of freedom of the microscopic 
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Lagrangian describing the complete system. The origin of the dissipation term and its quantum field theoretical 
treatment has been extensively discussed in the literature (for a recent review, see e.g. Ref. [131 )• However, as 
radiation is produced during the infiaton's evolution, the full treatment of the spectrum of perturbations no longer 
involves only that of the infiaton's perturbations but must also account for radiation perturbations. This makes 
the treatment of density perturbations in warm inflation similar to a multifiuid system. Since the larger is the 
dissipation the larger is expected to be the rate of radiation production, it becomes important the study of how 
the produced radiation and its perturbations backreacts on the infiaton's evolution and respective perturbations. 
In Ref. it was shown that as a consequence of this backreaction, the infaton's perturbations can grow as the 
dissipation of the infiaton increases. This increasing of the infiaton's perturbations with increasing dissipation 
can, therefore, severely constrain the model parameters in warm inflation so as to cope with the known measured 
results for the CMB radiation. This backreaction of the produced radiation on the infiaton's perturbations is 
larger the stronger is the coupling between the radiation and the infiaton; in particular the larger the power of 
T in the dissipative coefficient, the stronger the T dependence on the perturbations. This is exemplified by the 
results shown in Fig. [3] 

In this paper we have studied how this backreaction of the produced radiation, that can lead to this growth 
mode in the infiaton's perturbations, can be counterbalanced by the dissipative effects within the radiation 
fluid. Dissipative effects in the radiation fluid itself are described by viscosity terms. This is expected when 
the radiation fluid departs from equilibrium, which is the case in any dissipative system, where the produced 
radiation from the system does not immediately equilibrate in the radiation bath and its approach to equilibrium 
is controlled by viscosity coefficients, like the shear viscosity, the bulk viscosity and heat transport coefficients. 
We have focused on the dissipation effect as coming dominantly from a shear viscosity term in the ffuctuation 
equations. We have then shown that the shear viscosity can effectively damp the radiation fluctuations so as 
to avoid altogether the appearance of the growth mode in the resulting perturbations. The results we have 
obtained are model independent and we have shown that the overall effect of compensation of the growth mode 
and its control is determined by the ratio (sH/ {3pr), where Cs is the shear viscosity coefficient, H is the Hubble 
parameter and pr the radiation energy density. When the ratio is 0(1) or larger, the growth mode disappears 
completely. 

In this work we have only considered the coupling between the ffuctuations in the inflaton field with those in the 
radiation through the (temperature dependence on the) dissipation coefficient in the infiaton's dynamics. But in 
a thermal bath, the parameters of the infiaton's potential can also adquire temperature corrections. Even though 
these thermal corrections can be kept under control and small in Supersymmetry model building realizations for 
warm infiation [ij, E^, 113 , they can still be large enough to provide extra sources of couplings between infiaton 
and radiation fiuctuations and it should be interesting to analyze their effects in a future work. Likewise, there 
can be additional sources of dissipation in the radiation fluid, for example as coming from bulk viscosity, that 
can further help to damp any leftover growing modes as resulting from these additional couplings. In this work 
we have neglected the effects of the bulk viscosity on the grounds that it is in general much smaller than the 
shear viscosity. For example, for the self-interacting scalar field radiation discussed in section IV, the ratio of 
the bulk viscosity, C,b, with the shear viscosity for a high temperature radiation fiuid is Cft/Cs ~ thus, it 

is negligible for a weakly interacting radiation bath. But there may be other interactions and energy regimes for 
the radiation fiuid in which the bulk can be sizable and lead to effects in the density perturbation evolution (see 
e.g. Ref. [411]). 

All these effects, starting with the shear viscosity, will also impact the second order evolution of the pertur- 
bations and thus the calculation of the non-gaussinity. Forthcoming cosmological data are expected very soon 
to set the level of non-gaussinity of the primordial spectrum, which clearly will help to discriminate among 
infiationary models. Warm infiation, being a type of multi-fiuid model, falls into the category of models with a 
non-negligible value of the non-linearity parameter Jnl for non-gaussinity. This parameter has been computed 
for a T independent dissipative coefficient ^] , and recently the T dependence of T has been included [i^ , which 
provides an extra non- linear source in the field second order equation. However, if the coupling between field 
and radiation perturbations at first order is suppressed by viscous effects, qualitatively we expect the same to 
happen at second order. The question then is whether one simply recover the prediction for a constant dissipative 
coefficient, or non-linearities are further suppressed by viscous effects. These and other effects mentioned above 
will be studied elsewhere. 
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